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The fluorescence spectrum of a strongly pumped two level system is characterized by the Mollow 
triplet that has been observed in a variety of systems, ranging from atoms to quantum dots and su¬ 
perconducting qubits. We theoretically studied the fluorescence of a strongly pumped intersubband 
transition in a quantum well. Our results show that the many-electron nature of such a system 
leads to a modification of the usual Mollow theory. In particular, the intensity of the central peak 
in the fluorescence spectrum becomes a function of the electron coherence, allowing access to the 
coherence time of a two dimensional electron gas through a fluorescence intensity measurement. 


The fluorescence spectrum of a resonantly pumped two 
level system (TLS) is characterised by the three peak pro¬ 
file named after Mollow [1]. Under strong pumping, the 
levels of the TLS are split by the ac Stark effect into 
doublets, |±), whose splitting, O, is directly proportional 
to the pump amplitude. Of the four possible sponta¬ 
neous emission channels between two doublets, two are 
resonant with the bare frequency of the TLS, Co’ 12 , and 
the other two are at frequency CJ 12 ± U [see Fig. 1 for a 
scheme of the emission channels and frequencies]. From 
this simple picture it can be inferred that the emission 
from the central peak and from the satellites are in ratio 
1:2:1, a result that still holds for more refined theoretical 
approaches [2-4] . To date the Mollow triplet has been ob¬ 
served in a wealth of different systems well modelled by a 
TLS: atoms [5, 6], quantum dots [7-12], single molecules 
[13], superconducting qubits [14-16], and semiconductor 
quantum well (QW) excitons [17, 18]. 

Intersubband transitions (ISBTs) occur between two 
QW conduction subbands, with the lower one containing 
a two dimensional electron gas (2DEG) created through 
doping [19], temperature [20], or by optically exciting 
electrons from the valence band [21, 22]. Unbound ex¬ 
citations, ISBTs have a narrow absorption line thanks 
to the fact that conduction subbands are quasi-parallel 
[23, 24] (see Fig. 2 (a) for a graphical representation). 
Although in ISBTs the Mollow triplet has not yet been 
directly measured, the ac Stark effect, which can be in¬ 
terpreted as indirect evidence, has been clearly observed 
[25]. While ISBTs are usually theoretically described as 
collections of independent TLSs [25-28], it has recently 
been shown [29] that such an approximation breaks down 
in the nonlinear regime, when a macroscopic fraction of 
the total number of electrons is in the excited subband. 
As explained in Ref. [29], the origin of such a difference 
can be intuitively traced to the different dimensionality 
of the Hilbert space: 2^ for a collection of N TLSs and 

‘^N ) ISBT involving N electrons and neglect- 

. iv 

ing border effects, which tends to for large N. 

In a setup adapted to measure the resonance fluores¬ 
cence of the system, almost all the electrons in the first 
subband are excited at the same time [25] . We can thus 
expect that the TLS approximation, used to derive the 
theory of the Mollow triplet, will break down. In this 
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FIG. 1 . (a) 

strongly illuminated TLS giving rise to the Mollow triplet. 

(b) Emission lines from the dressed states occur at the bare 
frequency Ct;i2 and at the shifted frequencies Ct;i2 ± U. 


paper we will develop a predictive many-electron the¬ 
ory describing resonant fluorescence emission from ISBTs 
and other similar systems beyond the TLS approximation 
[29, 30]. We show that the physical picture is rather dif¬ 
ferent from the TLS case and that, while the three peaked 
structure is conserved, the relative intensity of the peaks 
becomes a function of the coherence of the electron gas. 
Our work thus hints at the intriguing possibility of mea¬ 
suring the coherence time of a 2DEG through a measure 
of its resonance fluorescence. In the following we intro¬ 
duce the theoretical description of a resonantly pumped 
ISBT and develop a perturbative theory of fluorescence 
emission that will allow us to calculate the intensity of the 
emitted radiation. We conclude by considering the lim¬ 
itations of our theory and the experimental constraints 
required to observe such an effect. The interested reader 
will find explicit analytic derivations in the Appendix. 

Eollowing the theory developed in Ref. [31], we con¬ 
sider two conduction subbands, the lower one containing 
a 2DEG, pumped by a coherent source resonant at the 
ISBT transition frequency between the two subbands, 
Co ’12 [see Appendix A]. Being the in-plane wave vector 
conserved in a planar structure, the pump, which has a 
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FIG. 2. (a) Energy dispersion of the first two conduction 

subbands in a doped QW as a function of the in-plane wave 
vector k. The first conduction subband is partially filled with 
a 2DEG up to the Eermi wave vector kp- (b) The ratio be¬ 
tween the emitted photon wave vector, q, and the electron 
wave vector uncertainty A/c, determines how many hnal states 
are accessible for a photonic emission. 


well-defined wave vector q, couples each electronic state 
of the first subband, with wave vector k, to a single state 
in the second one, with wave vector k + q, and vice versa. 
The set of all the electron levels in the two subbands is 
divided into coupled pairs, and the Hilbert space of the 
electronic system is consequently partitioned into four¬ 
dimensional subspaces. If only one of each pair of levels 
is occupied by an electron, then we can consider only two 
states in each of these subspaces, and the system exactly 
maps onto a set of pumped TLSs, each one having eigen¬ 
values huj±^\^ = in the referential frame rotating 

at the pump frequency, where Q is the Rabi frequency 
proportional to the pump amplitude. The relative eigen¬ 
states are 


l±> = ^ (^U+q ± ^Ik) |0el) ’ (1) 

where |0ei) describes the empty subbands and cj ^ is the 
creation operator for an electron in the jth subband, 
j = {1,2}, with in-plane wave vector k. As all interac¬ 
tions are spin-conserving, all sums over electronic wave 
vectors run implicitly also over spin. In order to span the 
whole Hilbert space of the system, we also have to include 
the states in which neither or both levels of the pair are 
occupied, |F, k) = 4,k+q4,k l^ei) and |^,k) = |0ei). As 
these states are either empty or Pauli blocked, they do 
not couple with the pump, and have eigenvalue zero in 
the rotating frame [31, 32], as shown in Fig. 3 (a). To 
each k value we can thus associate the four-dimensional 
Hilbert space: Hk = {|-, k), |F, k), |F;, k), |+, k)}. Af¬ 
ter having solved the Hamiltonian for the closed system, 
in order to describe its fluorescence, we need to couple 
it to the free electromagnetic field. Calling the cre¬ 
ation operator for a photon of energy huOq^q^, in-plane 
wave vector q, and out-of-plane wave vector the in¬ 
teraction takes the form 

^ ~ + H.C., (2) 

k,q,<7^ 


where Xq,qz light-matter coupling term, directly 

proportional to the ISBT dipole moment. As only trans¬ 
verse magnetic modes couple to ISBTs, we neglect pho¬ 
ton polarization [19]. Contrary to the TLS case, the 
interaction in Eq. (2) causes scattering between two- 
electron eigenvectors, i.e., involving in general states be¬ 
longing to HkG)Hk+q- These processes, initially studied 
in Ref. [31], are shown in Fig. 3 (b), for the two-electron 
initial state in Fig. 3 (c). 

There are four different scattering channels, two emit¬ 
ting at frequency uji 2 and the other two at While 

this could seem to lead to a result analogous to that ob¬ 
tained with standard TLS Mollow theory, there is a fun¬ 
damental difference: here the final state of all the emis¬ 
sion channels is the same, which is the product of a fully 
occupied and an empty state, |F, k) G) |F^, k + q). This 
opens up the possibility of having interference for the 
emission of the central peak, with its intensity depend¬ 
ing on the relative phases between the two electrons in¬ 
volved. A simple understanding of the phase-dependence 
of the emission intensity can be gained by considering the 
semi-classical picture of the two electrons cycling between 
the two subbands, driven by the pump. The interaction 
Hamiltonian in Eq. (2) leads to a diagonal process, in 
which one of the electrons scatters from its own upper 
state to the other electron’s lower state. This transition 
is Pauli blocked unless both electrons are at the same 
time in their upper state. The emission intensity will thus 
depend on the relative phase of the two cycling electrons. 

In order to make quantitative such a handwaving dis¬ 
cussion, we need to describe the full many-body state of 
the system undergoing the Rabi oscillations, taking care 
of the relative phases of the electrons, 0k, which will have 
a paramount role in the following. Under the approxima¬ 
tion that all the electrons participate [25], we can write 
this state as [see Appendix B] 


l^elW) = (g)^- L -^ (3) 


V2 

2,k+q 


= +cos(^J^)cJ J |0ei) • 


This state is not stationary and as such it cannot be used 
in the standard formulation of the Fermi golden rule to 
calculate the photon emission. An equivalent formulation 
can be obtained by calculating, to the first order in U, the 
total number of emitted photons per unit time, leading 
to [33], as shown in Appendix D, 

^ k,q,g. 

with /k,q(t) = - nt) - cos( ^’-~^‘-+y . In 

order to calculate the sums in Eq. (D4) we neglect spe¬ 
cific correlations between the phases, considering them 
as independent and identically distributed random vari¬ 
ables. This is a mean-field approximation that allows ns 
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FIG. 3. (a) Eigenvectors and eigenenergies of one of the spaces for one of the subspaces 7i\^. (b) Transitions in intersubband 

resonance fluorescence. At Ct;i 2 , two different scattering processes interfere and the intensity of the emitted light is given by 
their relative phases, 4cos^ ^k+q ^ Two-electron state in which the electrons oscillate between the two subbands with 
phases and 0k+q- 


to replace the sum over k with an average over the phase 
distribution, multiplied by the total number of electrons 

N. While an exact investigation of the phase dynamics 
is beyond the scope of the present work, we can give a 
qualitative description of phase diffusion. In a typical ex¬ 
periment the phases will be initially all equal, as all the 
electrons lie in the lower conduction subband. Once the 
continuous drive of the pump is turned on, the electrons 
will eventually diffuse with a coherence time Tcoh, leading 
to the experimentally observed wash-out of the Rabi os¬ 
cillations [34-36], see Appendix C for details. Under this 
approximation, and for times long enough to be able to 
resolve the peaks {t > ^), we can calculate the emission 
rate as [see Appendix D] 

r = I (1 + (cos(/))^ + (sin</>)^) ro(wi 2 ) 

+ i^ro(wi2 + O) + j^ro(a;i2 — fi), (5) 

where ro(ci;) is the standard dipole emission at frequency 
uj and the angular bracket is the average over the phase 
(j). While electrons are fully coherent, the phase distribu¬ 
tion will be strongly peaked around a well defined value 

O, and in Eq. (5) we obtain (cos0)^ + (sin0)^ = 1, inde¬ 
pendently of the specific value of 0. If we approximate 
ro(co’i 2 ±U) ro(co’i 2 ), the intensity of the central peak is 
in a ratio 1:4:1 with respect to the satellites. For t ^ Tcoh 
instead, the electron coherence is lost. The phases can 


in this case be considered as uniformly distributed and 
we have (cos0) = (sin0) = 0, recovering the usual ratio 
1:2:1 for the area of the three peaks and the expected 
asymptotic value for the total incoherent emission rate 
Fin = ^ [31] [see Appendix D]. 

A clear way to present this result is in terms of the 
many-electron Rabi oscillations. The fraction of electrons 
in the second subband can be calculated as [see Appendix 
C] 

n 2 {t) = (V’el(i)l IV’el(i)) = (sin 2 (^)), ( 6 ) 

which have a visibility decreasing exponentially in time 
due to decoherence (see the inset of Fig. 4), given by an 
envelope function normalized between 0 and 1, 

C = 2maxn2(^) — 1 = a/ (cos 0) ^ + (sin 0) ^. (7) 

Remarkably, Eq. (5) can be recast as a simple function 
of C, 

r = 1(1 + + i^ro(wi2 + n) + ^To{coi 2 - n),{ 8 ) 

which is plotted in Fig. 4, assuming that Fo(cji 2 ± U) ^ 

To (<^ 12 ), valid for co ’12 ^ U. From Eq. (8) we can clearly 
see how the relative intensity of the central peak de¬ 
creases from 4 to 2 with dephasing. 
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FIG. 4. The integrated intensity of fluorescence, plotted nor¬ 
malized to the fully incoherent value Fin, decreases as the vis¬ 
ibility C goes to zero. Inset: The collective Rabi oscillations 
of the second subband, n 2 {t) (black solid line, left axis) are 
washed out in a time TcoH- The visibility C (red dashed line, 
right axis) decreases accordingly in time. Here QtcoH = 20. 


Notice that, while the coherence time can be measured 
by four-wave mixing [37], or by the populations, through 
Eq. (6) [36], this requires one to perform either a pump 
and probe measurement, or to have an appositely crafted 
system with a shelving level [25]. Through F instead, the 
same information can be acquired via a time-resolved flu¬ 
orescence measurement. A main experimental challenge 
will be to discriminate the weak fluorescence signal from 
the pump. While polarization cannot be used to this 
aim, due to ISBTs selection rules, it is possible to exploit 
the broad angular distribution of the fluorescence. Sim¬ 
ilarly to what was done in Ref. [7], the sample can be 
engineered for the pump beam to be confined in a waveg- 
uided mode, while part of the fluorescence is emitted in 
non-guided modes, and collected outside of the sample. 

Above we implicitly assumed that the visibility C can 
be considered almost constant over a few Rabi oscilla¬ 
tions, in order to neglect the time dependency of Eq. (7) 
in Eq. (D4). Such a condition, which can be expressed 
as flTcoh ^ 1, can be fulfilled by present-day technol¬ 
ogy. The main contributions to Tcoh generally come from 
interface roughness scattering (IRS) and longitudinal op¬ 
tical phonon emission, which are both suppressed in wide 
QWs [19, 38-41]. Since IRS scales as the inverse sixth 
power of the QW length and the phonon emission is for¬ 
bidden when the intersubband gap becomes smaller than 
the optical phonon frequency [42] , coherence times of the 
order of a picosecond are achievable in wide QWs. In 
particular, in GaAs/AlGaAs QWs similar to the ones in¬ 
vestigated in Ref. [43], of bare frequency hwi 2 = 100 
meV, at T = 4.2 K, linewidths of 2.5 meV, correspond¬ 


ing to Tcoh = 0.5 ps, can be obtained even under intense 
pumping [25, 28, 44]. In such structures, Rabi splittings 
in excess of 10 meV are achievable for internal pump in¬ 
tensities of the order of ^MW/cm^ [25, 40], leading to 
flTcoh — 10. In Ge/SiGe QWs, the non-polarity of the Ge 
lattice inhibits Erohlich-mediated phonon emission and 
removes the constraint of performing measurements at 
cryogenic temperatures. While no observation of Rabi 
oscillations has been made in these systems, recent re¬ 
sults with linewidths of a few meV have been obtained 
up to T = 300 K for the bare frequencies, so Tcoh = 1 ps 
or longer is reachable by current technology [45]. Eor a 
structure with huJi 2 = 50 meV and Lqw = 10 nm as in 
Ref. [45], we expect flrcoh — 20. 

We developed our theory neglecting the Coulomb inter¬ 
action. While it is well known that in the case of parallel 
subbands its effects usually reduce to a renormalization of 
the intersubband transition energy [23] (the so-called de¬ 
polarization shift, vanishing for parabolic wells [46, 47]), 
Coulomb interaction could have a non-negligible impact 
in our case, as the presence of collective excitations [48- 
52] could spoil the independent electron picture we used 
[53] . Eor this reason we can consider our analysis as rig¬ 
orous only in the limit in which the depolarization shift 
is much smaller than the Rabi frequency, and plasmonic 
effects can be ignored. In both structures from Refs. [43] 
and [45] we estimated a depolarization shift of less than 
0.8 meV for a carrier concentration of n2DEG = 10^^ 
cm“^ [19, 21, 54, 55]. Our hypothesis of a depolariza¬ 
tion shift much smaller than the Rabi frequency is thus 
fulfilled for all but the most heavily doped structures. 

A last point worth stressing is that we considered a 
planar infinite 2DEG. While this is usually a very good 
approximation for high mobility samples, we still have 
to consider that the oscillating electrons are localized in 
the laser spot. If we call A/c the electron wave vector 
uncertainty due to such a confinement, the observation of 
interference of multi-electron scattering will require that 
the in-plane momentum of the emitted photon, g, obeys 
the relation q ^ Ak [see Eig. 2 (b) for a scheme of the 
emission process]. Eor a typical waist re 50 /im [25], 
this condition is well fulfilled already for huji 2 > 7 meV, 
but attention should be paid when applying our theory 
to narrow waists and THz QWs. 

In conclusion, we have shown that when ISBTs are 
strongly driven at resonance by an optical pump, the 
fluorescence has peculiar features that are not found in 
that of a collection of non-interacting TLSs. Although 
the TLS approximation is useful to describe Rabi oscilla¬ 
tions, the richer dynamics introduced by incoherent emis¬ 
sion processes requires a more sophisticated approach. In 
ISBTs the ratio between the three resonance fluorescence 
peaks deviates from the Mollow triplet due to interference 
effects that enhance the emission from the central peak. 
Thanks to this effect, the coherence time of a 2DEG could 
be accessed directly from a measure of the time resolved 
fluorescence intensity. 
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Appendix A: Driven Hamiltonian 

In this Apendix we give a brief review of the theo¬ 
retical framework, introduced in Ref. [31], which we 
use to study a strongly pumped intersubband transition 
(ISBT). The Hamiltonian of the pumped electronic sys¬ 
tem can be written as 

H = 

+ ? E (4.k+qCi,ke-*‘^-* + H.C.) , (Al) 

k 

where cj ^ is the second quantization operator creating 
an electron in subband j = {1,2} with in-plane wave vec¬ 
tor k and energy hujj^k^ ^ is the Rabi frequency propor¬ 
tional to the pump field amplitude, and cjl is the pump 
frequency, whose in-plane wave vector is q. We neglect 
to explicitly mark the spin index as all the interactions 
are spin conserving. In the following all the sums over 
electronic wave vector, which run until the Eermi wave 
vector, will thus implicitly sum also over spin. Being 
the photonic wave vector much smaller than the typical 
electronic one, and given that the conduction subbands 
are to a good approximation parallel, we can consider 
<^12 — ^ 2 ,|k+q| ~ ^1,/c- Passing in the frame rotating 
at the pump frequency, in the resonant case cjl = <^ 12 , 
we can reduce the Hamiltonian to the time-independent 
form 

^' = ^E(4,k+qCi,k + H.C.). (A2) 

k 

Each value of k in the sum in Eq. (A2) spans an inde¬ 
pendent four-dimensional Hilbert subspace, whose eigen¬ 
vectors are 

1 =^’ ~ ^72 ^^Ik+q ^ ^l,k) l^el) 5 

|F, k) = l^el) : 

|E;,k) = |0ei), (A3) 

the first two with energy and the last two with zero 
energy in the rotating frame. In Eq. (A3) the ket |0ei) 
represents the empty conduction band. 

Appendix B: State of the driven system 

We assume that initially the pump is “off” and that 
the ground state of the system thus describes a two di¬ 
mensional electron gas in the lower subband 

l^el(0)) = (g)4^k|0el). (Bl) 

k 


The state in Eq. (Bl) can be rewritten in terms of the 
two eigenstates of H' in Eq. (A3) as 

IV’ei(O)) = (g) I+, k) - I-, k)) , (B2) 

k 

where we added the relative phase of each oscillating elec¬ 
tron 0k- Those phases are all zero at t = 0, in order to 
recover Eq. (Bl), but they can evolve in time due to de¬ 
phasing processes, which cannot be accounted for in the 
Hamiltonian dynamics we considered. When the opti¬ 
cal pump is “switched on” the electrons start oscillating 
under the action of H'. The state in Eq. (B2) can be 
evolved using Eq. (A2) and Eq. (A3) as 

= (g)e*^[isin(^5^)4 k+q + cos(^^^)cJ ^] |0ei), 

k 

(B3) 

where in the second line of Eq. (B3) we expressed the 
state in terms of the fermionic operators acting on the 
empty subbands, which makes it evident that the elec¬ 
trons are oscillating between the two subbands. 

Appendix C: Visibility of the collective Rabi 
oscillations 

The fraction of electrons in the second subband can be 
calculated from the evolved electronic state in Eq. (B3) 
as 

'n2{t) = (t/’elWI E4.kC2.k 

k 

k 

(Cl) 

where N is the total number of electrons. In the last 
passage of Eq. (Cl) we assumed that the phases 0k are 
independent identically distributed (i.i.d.) random vari¬ 
ables, so that ^ ^(0k) = (^(0)) foi* any function ^(0). 

Assuming that Rabi oscillations are faster than dephas¬ 
ing, their visibility can be operationally quantified as the 
difference between the maximum and the average of the 
population oscillations, that is 

C = 2 ^maxn2(t) — , (C2) 

where we chose the normalisation factor in order to have 
0 < C < 1. Inserting Eq. (Cl) in Eq. (C2) and develop¬ 
ing the sine we obtain 

C = 2max[(sin2 

= 2max[—((cos^ cos^ ^ + (sin^ sin^ ^ 

+ |(sin0) sinf2t)] + 1. (C3) 
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FIG. 5. An illustration of phase diffusion over time, (a) At t = 0, all electrons lie in the first subband and the state is well 
described by Eq. (Bl), that is Eq. (B2) with all (/)k = 0. (b) Eor t TcoH the electrons oscillate coherently, (c) Eor t ^ TcoH 
the electrons start dephasing, (d) At t ^ TcoH coherence between electrons is lost and Rabi oscillations are averaged out. 


The maximum over time can be found by taking the time 
derivative of the expression between the square brackets 
in Eq. (C3) and by imposing it equal to zero, that is 

(sin0) cosf^t + ^(sin^ — (cos^ sinf^t = 0,(C4) 
which is satisfied for 


Ot = arctan 


(sin (j)) 


(cos2 I) - (sin2 f) 


(C5) 


After some straightforward algebra Eq. (C3) can thus be 
rewritten in its final form as 


C = y/ (cos 0) ^ + (sin (j)) ^. 


(C6) 


Let us discuss qualitatively how the electron coherence 
of the system changes over time in a typical experiment 
of continuous pumping. At t = 0 all electrons lie in the 
first subband and we can assume that all 0k = 0, as 
shown graphically in Fig. 5(a). As the pump is switched 
on, the electrons start oscillating and for t <C Tcoh rela¬ 
tive phases will not have diffused much, as shown in Fig. 
5(b). Their distribution will thus be strongly peaked 
around a well defined value 0, and the average values 
in Eq. (C6) will effectively reduce to evaluate the func¬ 
tions at 0: C = \/ cos^ 0 -h sin^ 0 = 1, for any value of 
0. At times comparable to the coherence time, t ^ Tcoh, 
with the pump still illuminating the system, more elec¬ 
trons will have been involved in incoherent scattering 
processes, as visible in Fig. 5(c), decreasing the visi¬ 
bility of the Rabi oscillations. Finally, for times much 
longer than the coherence time, t ^ Tcoh, the phases of 
the different electrons will be completely randomized, as 
schematized by Fig. 5(d). In this case the phase aver¬ 
ages in Eq. (C6) will be over an uniform distribution, 
(cos0) = (sin0) = 0, leading to (7 = 0. 


Appendix D: Photon emission rate 

In order to describe the fluorescence we need to intro¬ 
duce the electromagnetic field, whose free Hamiltonian 


can be generally written as 

-^ph = 

where is the creation operator of a photon with 

wave vector q in the plane of the quantum well and 
normal to it, whose energy is hujq^q^. The coupling be¬ 
tween the electromagnetic field and the ISBT is given in 
the rotating wave approximation by the interaction term 

^ ~ ^ H.C., (D2) 

k,q,g^ 

that in the rotating frame, and in interaction representa¬ 
tion with respect to iLph + H' ^ becomes 

k,q,g;z 

+ H.c. (D3) 


The rate of the emitted photons from the system, F, can 
be calculated by dividing the total number of emitted 
photons by the interaction time 

r = ’ (D4) 

ci,qz 


where |'0(t)), the total system wave vector in interaction 
representation, can be obtained by calculating to first 
order in V the evolution of the electronic state in Eq. 
(B2) and of the photonic vacuum, |0ph) 


IV’(^)) = IV’el(O))0|Oph) - T / F(t) |?/)el(O))0|Oph)dr. 

« J-t/2 


(D5) 


Using Eq. (D3) and Eq. (D5), the photon population 
can be calculated as 


Nph{t) = ^ JJ ^ dTdr' (t/’(0)| U(r) IV’(O)), 


(D6) 
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and after some straightforward algebra Eq. (D4) can be 
rewritten explicitly as 

r = E lx.,., / 


with /k,q(t) = cos( '^‘-+^‘-+^ - m) - cos( '^‘-y’-+-^ ). The 
integrand in Eq. (D7) contains only phase factors lin¬ 
ear in r and as such, for long enough times, Eq. (D7) 
will be given by a sum of squared Dirac delta-like func¬ 
tions. Using the usual trick commonly employed to 
derive the Eermi golden rule of formally transforming 
6{uj)‘^ = (5(0)(5(co’) = we obtain 

E lx.,.,P[4cos2(^Ji^|!s±a),5(wq,q^ - W12) 

k,q,g^ 

~ (‘^12 + ^)) + ^{^ q,qz ~ (“^12 ~ ^))]- 

(D8) 

As done previously for the electron visibility, we assume 
that the phases are i.i.d. random variables. The sum 


over k of the only term depending on k in Eq. (D8) thus 
becomes 

^E4cos"(^'^) 

k 

= ^ ^ 2 (1 + COS 0k cos 0k+q + sin 0k sin 0k+q) 

k 

^‘=‘ 2(1 + (cos0)^ + (sin0)^). (D9) 

Einally, inserting Eq. (D9) into Eq. (D8), and exploiting 
Eq. (C6), we obtain 

r = ^(1 + c^)^o{^i2) + i^ro(cji2 + u) + i^ro(cji2 — u), 

(DIO) 

with ro(co’) being the spontaneous emission rate of N two- 
level systems with bare frequency uj. Notice that, for long 
times t ^ Tcoh, C = 0 and we recover E 
as discussed in Ref. [31], the factor | comes from the 
fact that both the initial and the final state have average 
and uncorrelated filling factors 
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